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Abstract 

Based on the reported data until 18 March 2015 and numerical fitting via a simple formula of cumulative 
case number, we provide real-time estimation on basic reproduction number, inflection point, peak time and 
final outbreak size of ongoing Ebola outbreak in West Africa. From our simulation, we conclude that the first 
wave has passed its inflection point and predict that a second epidemic wave may appear in the near future. 
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1. Introduction 

Ebola virus disease (EVD) is a severe disease in 
humans which has infected nearly 25 thousand indi¬ 
viduals and claimed more than ten thousand deaths 
during the recent outbreak in West Africa, accord¬ 
ing to the report of World Health Organization dated 
18 March 2015 [1, 2]. The most affected countries 
are Guinea, Eiberia and Sierra Eeone. This study 
aims to provide some real-time estimations on the 
outbreak in these three countries using the reported 
cumulative case data. Specifically, we will estimate 
the following quantities: 

1. basic reproduction number 7?o, which is defined 
as the average new cases caused by a single in¬ 
fective individual during one infectious period; 

2. inflection point tc, which marks the time when 
the increment speed of cumulative case numbers 
starts to slow down; 

3. final outbreak size K, which indicates the total 
number of infectious cases throughout the out¬ 
break wave. 

4. peak time tp, which is defined as the critical time 
when daily infectious number reaches its maxi¬ 
mum. 
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All of these indicators provide quantitative informa¬ 
tion about severity of a disease outbreak. 


2. Methods 


Following [3], we study the epidemic model: 


^ s(t) + i(ty 
fismt) 


(1) 


i'(t) = 


s(t) + m 


ym. 


where S{t) and I{t) are the numbers of susceptible 
and infective individuals at time t, respectively. The 
constant jS denotes the transmission rate of the dis¬ 
ease, and the constant y corresponds to the removal 
rate of infective individuals. The basic reproduction 
number [4, 5] is given by 

Ro = -• ( 2 ) 

r 

It is noted that a disease outbreak occurs if and only 
if > 1. The differential system (1) can be solved 
explicitly and its solution is given by 

S(t) =i^[l +er(«o-l)(t-4)]-Ro/(Ro-l). 

I(t) =K[1 + er(«o-i)d-4)j-i/(Ro-i) (3^ 
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where K and tc are two eonstants of integration. Now, 
we define the eumulative infeetive ease number at 
time t as 


C{t) = 


r /3smt) 

loo S(t)+ /(O' 


(4) 


From (1) and (3), we have 


C(0 = K-K[1 + er(Ro-i)(f-(.)]-Ko/(Ro-i)_ (5) 


Here, the eonstant K = C(oo) has the biologieal 
meaning of final outbreak size. It ean be verified that 
C"(tc) = 0. Henee, tc is the infleetion point of C(t). 
We remark that the infleetion point tc is related to but 
different from another eommonly used quantity: the 
peak time, denoted by t^. The peak time is defined 
as the time when infeetive ease number aehieves its 
maximum, namely, /'(tp) = 0. It follows from (3) 


that 


tc + 


7 (^ 0 -!)■ 


( 6 ) 


In the ease when Rq is elose to 1, namely, ln7?o ~ 
- k we ean approximate the differenee tp - tc by 
1 /y. Thus, the peak time oeeurs about one infeetious 
period after the infleetion point [3]. 

Richards’ empirical model [6] was suggested to 
provide real-time estimation of a disease outbreak; 
see [7] for example. However, some of the parame¬ 
ters in Richards’ model do not have clear biological 
meanings [3]. The advantage of formula (5) is that 
all of the parameters in this formula have significant 
biological interpretations. We will use the explicit 
formula of C{t) in (5) to fit the reported cumulative 
case numbers of 2014 Ebola outbreak in West Africa 
and provide real-time estimation of basic reproduc¬ 
tion number Rqj inflection point tc, final outbreak size 
K and peak time tp. 

As pointed out in [3], one should fix the value of 
y, the removal rate of infective individuals, to re¬ 
solve possible overfitting problems. Note that 1/y 
can be regarded as the infectious period which char¬ 
acterizes the average duration of an individual be¬ 
ing infective. In most cases, an individual is re¬ 
moved from the infective group either by recovery 
or death. For the fatal cases of Ebola virus disease, 
death usually occurs between 6 and 16 days (with 
mean 7.5 days) after onset of symptom; and for the 
non-fatal cases, patients may improve their symp¬ 
toms at around day 6 but need more time to recover 


[8]. Convalescent patients may still be infective be¬ 
cause the Ebola virus RNA may remain in the body 
fluid for a couple of weeks even though the risk of 
transmission from them is low [9]. It is thus reason¬ 
able to assume the infectious period to be 7.5 days 
with some possible perturbations in the interval be¬ 
tween 6 and 16 days. In our simulation, we first fix 
1/y = 7.5 days to estimate the basic reproduction 
number Rq, inflection point tc, peak time tp and final 
outbreak size K using reported cumulative case data 
of Ebola virus in Guinea, Eiberia and Sierra Eeone, 
respectively [1, 2]. Also, we provide the 95% con¬ 
fidence intervals of each estimated parameter value 
using bootstrap method. Next, we vary the value of 
the parameter 1 /y from 6 to 16 days and investigate 
the sensitivity of fitted parameter values. 

3. Results 

The basic reproduction number Rq is estimated as 
1.116(95%CI: 1.115-1.116) for Guinea, 1.226(95% 
Cl: 1.225-1.228) for Eiberia, and 1.181 (95% Cl: 
1.181-1.182) for Sierra Eeone. The inflection point 
tc is estimated as 21 November 2014 for Guinea, 24 
October 2014 for Eiberia, and 28 November 2014 for 
Sierra Eeone. As shown in Table 1, the lengths of 
95% confidence intervals for the estimated inflection 
points are no more than one day. 



Figure 1; Fitted graph for the reported cumulative cases in 
Guinea. The dots are real data and the curve is plotted using 
fitted results. 

The fitted curves together with reported cumula¬ 
tive case data are illustrated in Figure 1 (Guinea), 
Figure 2 (Eiberia) and Figure 3 (Sierra Eeone). It 
is noted that in each of these three figures, there is 
a jump on the reported cumulative case numbers in 
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final outbreak size K 

basic reproduction number Rq 

inflection point tc 

peak time tp 

Guinea 
Liberia 
Sierra Leone 

3268 [3257, 3274] 
8630 [8605, 8660] 
11227 [11198, 11253] 

1.116 [1.115, 1.116] 

1.226 [1.225, 1.228] 

1.181 [1.181, 1.182] 

266 [265, 266] 
238 [237,238] 
273 [272, 273] 

273 [272,273] 
244 [244, 245] 
279 [279, 280] 


Table 1: Estimated parameter values with 95% confidence intervals. Here, day 1 corresponds to 1 March 2014. So, days 266, 238 
and 273 correspond to 21 November 2014, 24 October 2014, and 28 November 2014, respectively. 



Figure 2; Fitted graph for the reported cumulative cases in 
Fiberia. The dots are real data and the curve is plotted using 
fitted results. 



Figure 3; Fitted graph for the reported cumulative cases in 
Sierra Feone. The dots are real data and the curve is plotted 
using fitted results. 

late October 2014. This is due to a more compre¬ 
hensive assessment of patient databases on the World 
Health Organization report dated 29 October 2014 
[10]. Among these three countries, Liberia has the 
most significant gap, which may account for the re¬ 
sult that the inflection point for Liberia is about one 
month earlier than the other two countries. 

We also fit the final outbreak size as 3268 for 
Guinea, 8630 for Liberia, and 11227 for Sierra 
Leone. All of these estimated values are smaller 
than cumulative case numbers reported on 18 March 
2015. This indicates that another potential outbreak 
wave may be approaching [7]. 


Now, we regularly increase the value of infectious 
period l/y from 6 to 16 days, and conduct numerical 
simulations. It is noted that the fitted values of ba¬ 
sic reproduction number Rq will also increase from 
1.092 to 1.252 for Guinea, from 1.179 to 1.503 for 
Liberia, and from 1.144 to 1.400 for Sierra Leone; 
see Figure 4. The estimated final outbreak size stays 
in a range of [3261, 3307] for Guinea, [8620, 8677] 
for Liberia, and [11209,11319] for Sierra Leone; see 
Figure 7. On the other hand, the inflection point L 
and peak time tp do not vary too much; see Figures 5 
and 6. For Guinea, decreases from 266 (21 Novem¬ 
ber 2014) to 264 (19 November 2014). For Liberia, L 
decreases from 238 (24 October 2014) to 235 (21 Oc¬ 
tober 2014). For Sierra Leone, L decreases from 273 
(28 November 2014) to 270 (25 November 2014). 
The peak time tp increases from 272 (27 November 
2014) to 279 (4 December 2014) for Guinea, from 
244 (30 October 2014) to 248 (3 November 2014) for 
Liberia, and from 279 (4 December 2014) to 284 (9 
December 2014) for Sierra Leone. We observe that 
the fitted inflection point and peak time tp are sta¬ 
ble under perturbations on the infectious period l/y. 



Figure 4: Estimated values of basic reproduction number when 
infectious period increases from 6 to 16 days. 
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Figure 5: Estimated values of inflection point when infectious 
period increases from 6 to 16 days. Here, day 1 corresponds 
to 1 March 2014. So, days 266, 238 and 273 correspond to 
21 November 2014, 24 October 2014, and 28 November 2014, 
respectively. 



Figure 6: Estimated values of peak time when infectious pe¬ 
riod increases from 6 to 16 days. Here, day 1 corresponds to 
1 March 2014. So, days 266, 238 and 273 correspond to 21 
November 2014, 24 October 2014, and 28 November 2014, re¬ 
spectively. 

4. Discussion 

This study provides real-time estimation of basie 
reproduetion number, final outbreak size, inflection 
point and peak time for the onging Ebola outbreak in 
West Africa using reported cumulative case data. 

The fitted basic reproduction numbers are smaller 
than those estimated in [11] where only the data until 

20 August 2014 was used. This indicates that the 
disease control policy became more effective during 
the late stage of the outbreak. 

We also observe that the increment speed of cu¬ 
mulative case number began to slow down after 

21 November 2014 in Guinea, 24 October 2014 in 
Liberia, and 28 November 2014 in Sierra Leone. 
The estimated inflection points for Guinea and Sierra 
Leone are close to each other, but the one for Liberia 
is about one month earlier. This is due to a signif¬ 
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Figure 7: Estimated values of final outbreak size when infec 
tious period increases from 6 to 16 days. 


icant increase on reported cumulative case number 
dated 29 October 2014; see Ligure 2 and [10]. 

Lrom Table 1, we note that the estimated peak time 
has about one week’s delay after the estimated inflec¬ 
tion point, while the infectious period is fixed as 7.5 
days. This supports the conclusion in [3] that the 
peak time occurs about one infectious period after 
the inflection point. 

If we vary the infectious period from 6 to 16 days, 
the estimated basic reproduction number stays in a 
range of [1.092, 1.252] for Guinea, [1.179, 1.503] for 
Liberia, and [1.144, 1.400] for Sierra Leone. The es¬ 
timated final outbreak size ranges from 3261 to 3307 
for Guinea, from 8620 to 8677 for Liberia, and from 
11209 to 11319 for Sierra Leone. The estimated in¬ 
flection point and peak time are much stabler and 
only varies within a small interval. This demon¬ 
strates that our method has a significant accuracy in 
capturing the inflection point and peak time. 

The values of final outbreak sizes in three coun¬ 
tries are all underestimated, which can be considered 
as a warning signal of a second outbreak wave. 

Appendix A. Explicit solution of (1) 

Here, we provide the detail in solving the system 
(1). Lirst, we add the two equations in (1) to obtain 
S'{t) -I- /'(O = -y7(0- Coupling this with the first 
equation of (1) yields 

d{S + I) _ y{S + 7) 
dS “ fiS ■ 

This equation is separable and its solution can be 
written as 

Cl (5 +7) = 


4 

















where Ci > 0 is a eonstant of integration. Next, 
we use the above relation to eliminate I in the first 
equtaion of (1). It follows that 

S'{t) = 

This is a Bernoulli equation. We set u(t) = [S 
and obtain 

u'(t) = (/3-y)[u(t)-Ci]. 

The above equation becomes separable and its solu¬ 
tion is given by 

u(t) = Ci[l + C 2 e^-^^*] 

with C 2 > 0 being another constant of integration. 
Now, we have 

S(t) = 

Substituting this into the relation between S and I 
gives 

i(t) =cf ^’'“^^[1 + 

Finally, we set K = and tc = -(In C 2 )lifi - 7 ) 

to rewrite the formulas of S (t) and I(t) as 

S(t) =K[l + e(^-r)(»-f.)]/3/(r-/J); 

I(t) =K[l + e(^-r)(»-f.)]r/(y-/3) 

_ ^[1 -I- g(/3-r)(f-fc)j/3/(y-/3) 

This is equivalent with (3) in view of (2). 
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